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ABSTRACT 

The Cataclysmic Variable (CV) population harbors a diverse range of donor stars and accreting 
white dwarfs (WDs). A range of WD masses is expected, from low mass Helium core WDs, to massive 
WDs which have previously accreted at rates high enough for Hydrogen to burn steadily. Furthermore, 
a wide range of Helium enrichment is expected in the accreted material depending on the degree to 
which the donor star is evolved. We investigate the impact of this diversity on the range of effective 
temperatures (Toff) for which g-modes are unstable. Motivated by earlier theoretical studies, we use 
a simple criterion for g-mode excitation: that the thermal time at the base of the convection zone 
becomes longer than (some multiple of) a fiducial shortest g-mode period. The critical Toff below 
which modes are unstable ("blue edge") then depends on both surface gravity, g, and He abundance, 
Y . The Hydrogen/first Helium ionization instability strip is more sensitive to g than Y . We find that 
(for solar composition envelopes), relative to a fiducial WD mass 0.6 Af©, the blue edge for a 0.4 M© 
He core WD shifts downward by « 1000 K, while that for a massive w 1.2 M© WD shifts upward by 
« 2000 K. The second Helium ionization instability strip exhibits strong dependences on both g and 
Y. Surprisingly, increasing Y by only 10% relative to solar creates an instability strip near 15,000 K. 
Hence CV's below the period gap with evolved donor stars of F > 0.4 may have an "intermediate" 
instability strip well outside of the isolated DA and DB variables. This "intermediate" instability 
strip also occurs for low mass He WD with solar composition envelopes. The lack of pulsations in 
CV's with Teff in the pure Hydrogen ZZ Ceti instability strip is also easily explained. 
Subject headings: binaries: close — gravitational waves-novae, cataclysmic — variables — white dwarfs 



1. INTRODUCTION 

As they cool in isolation, WDs cross regions of pulsa- 
tional instability and undergo non-radial oscillations, in 
particular g-modes (see Gautschv & Saio 199q[). Those 
with pure H atmospheres (DAV) cross the ZZ Ceti insta- 
bility strip when Tcff « 11, 000 - 12, OOOK (see Mukadam 
et al 2004; Gianninas et al. 2005), whereas those with 
pure He envelopes ( DBV) become unstable when Teg « 
22,400 - 27,800 K l|BeauchamD et al.llT999|l . The com- 
position of these envelopes is rather pure, and the range 
of WD masses is that expected from isolated stellar evo- 
lution, favoring M w Q.QMq. Careful analysis of the 
adiabatic pulsations for these WDs yield accurate mea- 
suremen ts of the WD masses (M), shell m asses, and spin 
(see e.g. IBra,d1evll2?i?Ill IKenler et al .11200^ . 

The recent discoveries of similar period pulsations from 
accreting WDs in CVs provides a window to their inter- 
nal properties. The accretion of > O.IM© over 10^ years 
impacts both the WD's rotation rate and mass, param- 
eters easily diagnosed through seismic studies. Starting 
with the initial discovery (jvan Zvl ct al. 2000) of 3 os- 
cillation peri ods in the WD primary of GW Lib, there 
are now 10 IWarner fc WoudtJ IpM IWondt fc Warner 
20041: lArauio-Betancor et al.1 l2005t iPatterson et al 



22051; iSzkodv et al.' '2005^ 'Vanland ingham et all l200'5r 
Qansicke et al. 2006) pulsating WDs in CVs. These are 
all found in CVs below the period gap, where a quiescent 
spectrum shows e vidence for a WD with Tes < 25, 000 K 
l|Sionlll9911 I1999D Calculations of WD hea ting by pro- 
longed accretion IjTownslev fc Bildstenll200^ explain the 



Toff's at these orbital periods, and were employed for the 
first adiabatic analysis of GW Lib, yielding both the WD 
and accreted layer mass (Townslcv ct al. 2004). 

We show here that the expected diversity of WD 
masses and accreted layer compositions (in particular, 
enriched He abundance from an evolved donor) in these 
systems will lead to a broad range of Toff's where WDs 
will pulsate. For example, low-mass He core WDs can 
pulsate at T^s ~ 10, 000 K, well below the ZZ Ceti insta- 
bility strip. The enhanced Helium abundance expected 
from an evolved secondary can also push the instability 
strip to values as high as 20, 000 K. These effects likely 
explain why many of the observed pulsators in CVs are 
"outside" of the conventional Z Z Ceti in stability strip 
appropriate to pure H envelopes IjSzkodv et al...2002,) . 

In §2, we clarify the influence of WD mass and envelope 
composition on mode driving. Assuming that Brickhill's 
convective driving mechanism operates, we explore the 
depth of the convection zone as a function of composition 
and surface gravity. A near-solar mix of elements allows 
the possibility of both H/Hel and Hell ionization zones, 
opening up additional frequency ranges for instability of 
hotter and/or rapidly rotating WD. Section 3 contains a 
discussion of the expectations for the occurrence of He 
core WDs in CVs, and explains why many CVs can have 
an evolved donor that provides He rich material to the 
WD envelope. We close in §4 by summarizing our new 
understanding and discussing future work. 
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Accreting WDs consist of a geometrically thin, mostly 
non-degenerate accreted envelope overlying a degener- 
ate core. The envelope composition (of mean molecular 
weight /i) is inherited from the donor star, while the core 
composition can be Helium for low mass (< 0.45 Mq) 
WDs or a mixture of Carbon and Oxygen for larger 
masses (0.6 — I.IMq). A luminosity from "compressional 
heating" of accretion is generated deep in the envelope 
and core (Townsley & Bildsten 2004), hence the luminos- 
ity is nearly constant at the shallow depths where mode 
driving (and damping) occurs. 

The composition discontinuity at the base of the ac- 
creted envelope (where the temperature is Ti,) sepa- 
rates the WD into two resonant cavities: the envelope 
and the core (Townsley et al. 2004). The highest fre- 
quency gravity wave (g-mode) lives mainly in the enve- 
lope (and has no radial nodes there) and has a frequency 
w ~ [gHf/^kx - 27r/(100 s) (for I = 2), where H is 
the scale height at the base, (gHy^^ ~ {kBTi,/iim,py^^ 
is roughly the sound speed at the base, and k± — 
[l{l + 1)]^'^/R is the horizontal wavenumber. Internal 
gravity waves trapped in the envelope have frequencies 
which are smaller by w 1/n, where n is the number of 
radial nodes there. 

Gravity waves are driven overstable by converting part 
of the escaping heat flux into mechanical motion. This 
occurs near ionization zones, either due to a rapid out- 
ward increase in opacity in radiative zones ( "the k mech- 
anism" , Dziembowski 1977) , or more likely, due to a 
convection zone caused by such rapid increase in opac- 
ity (Brickhill's "convective dr iving" mechanism. iBrickhill 
IT9831 IWu fc GoldreichllT99l . We focus on BrickhiU's 
mechanism, where the key quantity is wTth.bcz, the prod- 
uct of the mode frequency and the thermal time, Tth,bcz, 
at the base of the conv ection zone. When wrthjbcz ^ 1/4 
IjWu fc Goldreichlll999|) . the thermal adjustment time in 
the convection zone is longer than the mode period, and 
flux perturbations are absorbed by the convection zone, 
acting to put energy into the oscillation. In this pa- 
per, we adopt the instability criterion P < 87rrth,bcz, 
where P = 100 s is the fiducial shortest period g-mode. 
We will not calculate the "red edge" Tcff of the insta- 
bility strip, below which the pulsations arc unobserv- 
able. For isolated DAVs, the i nstability strip is quite 
narrow (« 1000 K; iMukadam et al. 2004 j, and is much 
broad er (w 5000 K) for the DBVs (,BeauchamD et all 
11999(1 . Lastly, when two convection zones are present, we 
only consider driving by the inner convection zone as this 
most closely resembles the theory deriv ed for the DAVs 
IBrickh ill IQSa iWu & GoldreichI 119991) . However, driv- 
ing may also be possible by the outer convection zone, 
if the radiative layer between the two convection zones 
does not contribute too much damping. We leave this 
issue for future studies. 

The variation of Tth,bcz with g and Y is found by con- 
structing atmosphere models that use the OP AL opacity 
and equation of state ijlglesias fc R,ogerJl99(itl with solar 
metallicity {Z = 0.02). These are plane parallel, constant 
gravity and flux envelopes, using the ML2 mixing length 
prescription ( Bergeron et alJ ll992) in convective regions, 
in which the mixing length is set equal to the pressure 
scale height. 

Figure ^ shows the longest driven mode period. 



87rrth,bcz, as a function of T^s for g = lO'* cm s^^ and 
a range of y. At a given T^s there can be two distinct 
convection zones: an outer convection zone due to H 
and first He ionization (H/Hel) at short thermal times, 
and an inner convection zone due to second He ioniza- 
tion (Hell) at long thermal times. The six different lines 
show y = to F = 0.98, from left to right. For solar 
abundance {Y — 0.28, leftmost dashed line), we find that 
87rrth,bcz exceeds 100 s for T^s < 12, 000 K; the blue edge 
is not raised significantly by addition of a solar fraction 
of He relative to a pure H atmosphere. At solar He frac- 
tion, the inner Hell convection zone is only present at 
high temperatures (> 22,000 K) and only drives £ ^ 1 
modes with periods P < 10 s, which are unobservable. 
However, increasing Y above solar, even by a moderate 
amount, has a dramatic effect on the Hell convection 
zone. For Y = 0.38 (solid line), the inner Hell convec- 
tion zone is present down to Tcff ~ 14, 000 K, and raises 
the blue edge to near 16,000 K. We call this the "inter- 
mediate" instability strip. 

The Y = 0.38 line also indicates a feature of mode driv- 
ing in CV WDs, in which a single object can have two 
instability strips. The upper strip in Teff, arising from 
the Hell convection zone, has a red edge determined by 
the truncation of the Hell convection zone at lower Toff, 
rather than BrickhiU's (1983) red edge due to vanishing 
flux perturbation at the surface. Instability can then 
reappear at T^s S 12,000 K. This only occurs over a 
limited range of Y. At F = 0.48 (middle dashed line), 
the Hell convection zone is not truncated, and below the 
blue edge it extends across 87rTth — 100 s. Rather than 
the H/Hel convection zone rising up to create another in- 
stability strip, the two convection zones merge when the 
radiative region between them disappears. The presence 
of H will, however, create an extremely wide instability 
strip by keeping the convection zone thin compared to 
the Y = 0.98 case (right dot-dashed line). 

A similar experiment is shown in Figure |21 where g 
varies from g = 10^*' to 10*'^ cm s~^ keeping solar com- 
position (X = 0.7, Y = 0.28, Z = 0.02). A gravity of 
g = 10"^-^ cm s"^ corresponds to a M w 0.4AfQ WD, 
while g = 10* cm s~^ gives M w 0.6Mq, and g — 10^ cm 
s~^ gives M Ki 1.2Mq. As in Figure ^ the blue edge is 
the intersection of a curve with the horizontal line repre- 
senting the fiducial shortest period g-mode. As discussed 
by Arras et al. (200J|), the blue edge for the H/Hel ion- 
ization zone moves up (down) by « 2000 K for a factor 
of 10 increase (decrease) in gravity. Hence, a blue edge 
range of w 4000 K is expected from the lowest mass He- 
core WD to the massive > I.IMq Carbon/Oxygen WD. 
In short, the range of WD masses expected in CVs should 
extend the H/Hel instability strip to be both well above 
and below the isolated WD ZZ Ceti instability strip. For 
the Hell ionization zone, note the appearance of the "in- 
termediate" instability strip at « 15, 000 K for low mass 
WDs. Since Tth.bcz decreases as g increases, the blue edge 
for Hell ionization decreases with higher g. This behav- 
ior is opposite that found in the nearly pure He case, 
where Tth,bcz increases as g increases (see figure |31). 

To elucidate the sensitivity of the blue edge to the wide 
range of gravities and enrichments possible for CV WDs, 
Figure 01 shows contours of constant blue edge Tcs in 
logg and Y space. At small Y, driving occurs due to the 
H/Hel convection zone, while at large Y driving is due to 



Pulsational Instability in Accreting WDs 



loS 



10 E 



I lo' 

o 10 

•a 

Oh 
O 

H-1 

10 



10 



-1 



10 



1 — I — I — I — I — I — I — I — r3 




J I I I 



TffdO K) 



Fig. 1. — Longest period driven mode, S-TrTth ^czi vs. Tefj for different H to He ratios, but fixed gravity g = fO cm s~ and metallicity 
Z = 0.02. A fiducial sfiortest period g-mode, shown at 100 s, is driven when the curve is above the horizontal line. There can be two 
convection zones associated with second He ionization (Hell, upper branches), or first He ionization and/or H ionization (H/Hel, lower 
branches). Lines are shown for Y = 0.0 (left dot-dashed), solar composition (Y = 0.28, leftmost dashed), Y = 0.38 (solid), 0.48 (middle 
dashed), 0.58 (rightmost dashed), and 0.98 (right dot-dashed). Dotted lines only connect related curves and do not represent convection 
zone boundaries. 



the Hell convection zone. For low and moderate g, the 
transition occurs abruptly when the Hell convection zone 
extends above 100 s, causing the blue edge to jump up 
by several thousand K. The dashed line indicates the Y 
above which the Hell convection zone is always present. 
As discussed above, there is a small region, indicated on 
the plot, where there are two separate instability strips. 
Dotted lines show the blue edge of the lower instabil- 
ity strip. We find that starting from solar abundance 
{Y = 0.28, dot-dashed Hue), g must be increased signif- 
icantly, half an order of magnitude, in order to increase 
the blue edge T^s by 1000 K. In contrast, a modest in- 



crease in Y can effect a very sharp increase in the blue 
edge. There are thus two types of accreting pulsators, 
those with blue edge Tcff > 15, 000 K that are slightly 
enriched, and those closer to solar composition with blue 
edge Tcff < 15, 000 K, more similar to normal ZZ Ceti 
stars. Observationally, these two instability strips may 
overlap in Tgff depending on the location of the red edge. 

3. HELIUM ENRICHMENT IN EVOLVED CATACLYSMIC 
VARIABLES 

Both theory and observation (see Townsley & Bildsten 
2003) tell us that non-magnetic CVs need to be below 
the period gap (Porb < 2 hr) to reach the Toff < 25, OOOK 
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Fig. 2. — Longest period driven mode, STrrth bczi v! 
gravity for each line is, from left to right at low T^f{, 
and do not represent convection zone boundaries. 



■ Tcff for different gravity, but fixed (solar) composition Y = 0.28 and Z = 0.02. The 
g = lO'"'^, lO'"'^, 10^'^, and 10*'^ cm s"'^. Dotted lines only connect related curves 



range where pulsations are possible. This is indeed where 
the known pulsating WDs in CVs occur. We showed in 
the previous section that, due to the dependence on Y 
and g, the instability strips for accreting WDs will be 
richer than in isolated DA or DB WDs. 

The first diversity to consider is the WD mass. The 
expectation (see Howell, Nelson & Rappaport 2001; here- 
after HNR) is that WDs below the period gap will have a 
large range of masses, all the way down to low-mass He 
core WDs (M < 0.45Mq). HNR's population synthe- 
sis calculations showed that as many as 20% of the CVs 
with Porb < 2 hr will have low mass He WDs. These will 
have surface gravities much lower than any of the iso- 
lated DAVs, for which the recent tabulations (Mukadam 



et al 2004; Gianninis et al. 2005), show no DAV with 
log 5 < 7.7 - 7.8. The pulsator HS 2331+3905 h as a low 
Toff' ~ 10,500 K l|Arauio-Betancor et al.ll200l . which 
from Figure |21 is quite reasonably inside the instability 
strip for solar abundance with logg < 7.5. The inferred 
low mass of 0.4 Mq makes HS 2331+3905 a candidate 
He core WD. 

The instability strip is also sensitive to the He abun- 
dance in the WD envelope. At accretion rates > 10"^'* 
Mq yr^i for Toff > 11, 000 K, there is no t ime for 
relative H/He separation IjPaauette et al.lll986|) . Since 
even the small amount of quiescent disk accretion ex- 
ceeds this rate, the envelope He abundance is set by the 
donor star. If the donors were unevolved, then the range 
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Fig. 3. — Contours of constant T^s for the blue edge of the instabiUty strip in gravity, logg, and HeUurn mass fraction, Y. The blue 
edge is approximated by the Teg for which Sttth^ bcz = 100 s. For large Y, driving is due to the Hell convection zone, while for small Y it 
is due to the H/Hel convection zone. The dashed line indicates the Y above which the Hell convection zone is always present. The sudden 
increase of the blue edge with Y occurs where the HeH convection zone reaches above 100 s. Between this and the dashed line, there will 
be two separate instability strips caused by the absence of the HeH convection zone at some Tiff's. The blue edge of the lower instability 
strip is shown by the dotted lines. Solar composition, Y = 0.28, is indicated by the dash-dotted line. 



would be narrow, Y — 0.25 — 0.3, only depending on 
the initial nietallicity of the donor. However, Pylyser 
& Savonije (1989) raised the distinct possibility of initial 
donor masses for CVs with M > IMq, which, depending 
on the age of the system at the onset of mass transfer, 
can have undergone nuclear burning in the core to raise 
the He abundance. Such donors are required to explain 
supersoft sources as stably burning WDs (van den Heuvel 
et al. 2002). By the time such donors are below the pe- 
riod gap and fully convective, such a core He enrichment 
will be evident at the surface and in the accreting matter. 
This gives reason to explore the range ofY = 0.3 — 1.0. 
The prevalence of evolved donors amongst CVs with 



-Porb < 8 hours depends on many factors, including the 
efficiency of the common envelope phase, the outcome of 
the rapid thermal timescale mass transfer expected at the 
onset of Roche lobe filling, and the nature of magnetic 
braking in the pre-CV stage. (Schenker et al. 2002; Pod- 
siadlowski, Han & Rappaport 2003; Andronov & Pinson- 
neault 2004; Kolb & Willems 2005). Podsiadlowski, Han 
& Rappaport (2003) performed a binary population syn- 
thesis calculation that allowed for evolved donors with a 
range of initial He abundances amongst initially massive 
(M > Mq) donors. They found that these could actually 
be the dominant population for Porb > 5hours, but that 
the large number of CVs injected at shorter orbital pe- 
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riods dilutes their fraction at shorter periods. Below the 
period gap, between one in ten and one in three would 
have evolved companions with Y > 0.5 The previously 
discussed observables for evolved donors are a later spec- 
tral type (at a fixed orbital period), and exposed mate- 
rial that is Nitrogen enhanced and Carbon poor due to 
the CN cycle. This latter hypothesis has been tested by 
Gansicke et al. (2003) through measurements of UV line 
ratios in a diverse set of CVs, both above and below the 
period gap. They concluded that as many as 10-15% of 
their sample could have evolved donors. 

We are introducing here an additional probe of evolu- 
tion, which is the impact of the He abundance on the in- 
stability strip. More evolved donors will show pulsators 
at hotter T^s than those which are unevolved. In the 
absence of more detailed population study predictions of 
the He abundance (typically calculated, but not reported 
in depth since, until now, no observational probe of He 
abundance was available), we can say very little. Due to 
mass transfer stability considerations, the donor stars for 
a He core WD should be initially low mass and therefore 
unevolved. Hence the correct low gravity models to con- 
sider are those with Y = 0.25 — 0.3, but no higher. In the 
opposite sense, we would expect highly evolved donors to 
have preferentially more massive WDs. Whether these 
correlations lead to a range of Teff 's or other correlations 
remains to be seen. 

4. CONCLUSIONS AND FUTURE WORK 

We have shown that detecting pulsations in accreting 
WDs below the period gap has the promise to reveal the 
diversity expected in both the WD masses and the ac- 
creted He abundance. We do not expect one instability 
strip as known for isolated DA WDs, but rather a reveal- 
ing diversity of Tcff 's. In particular, hot CV primaries like 
that in GW Lib with Toff ~ 14, 000 K can pulsate due 
to enhanced Helium abundance, or extreme surface grav- 
ity (large or small). Cool CV primaries like that in HS 
2331-^3905 with Toff = 10, 500 K, can pulsate due to low 
surface gravity. Both of these possibilities are naturally 
present in the diversity expected within the population 
of CVs. Indeed, some systems can have two separate 
instability regions in T^s. 

Little is known with certainty about the WD mass dis- 
tribution, but evolutionary models provide a reasonable 
expectation for the level of Helium enrichment of the 



accreted material. CVs in which the companion is the 
remaining core of a slightly evolved star can have He sig- 
nificantly enriched above solar. Below the period gap, 
between one in ten and one in th ree CVs should have 
evolved co mpanions with Y > 0.5 IjSchenker et alJl2002t 
fPodsiadlowski et al., 2003') . 

Important work which remains is validation of our as- 
sumption that the numerical factor relating rth,bcz and 
the driven mode period does not vary significantly with 
either g or Y. A large variation (more than unity) is not 
expected, but quantifying this statement is necessary for 
using the presence or absence of modes to constrain prop- 
erties of particular systems. Additional damping due to 
the opacity effects of metals in the radiation zone may 
affect this prefactor, shifting the blue edges evaluated 
here to lower Toff. Inefficient convection in surface lay- 
ers would decrease the rate of convective driving, again 
causing the blue edge to shift downward. 

The WDs in CVs are expected to be rotating rapidly 
due to the angular momentum gain from long term disk 
accretion. This can qualitatively chang e the frequency 
spectrum of non-radial oscillations. Dzie mbowskil 1^97^ 
showed that g-modes with large numbers of angular 
nodes, and hence large frequency, can be driven even 
in relatively hot WDs by the K-mechanism acting in the 
He-II ionization zone. Such modes are unobservable in 
non-rotating stars due to averaging of the flux perturba- 
tion over the stellar disk. Rotationally modified g-modes, 
which are squeezed into an equatorial band, may be both 
unstable and observable, since they can have large hori- 
zontal wavenumber (and frequency) for a small number 
of angular nodes. 
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